Electric-octupole and pure-electric-quadrupole effects in soft-x-ray photoemission 
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Second-order [O(fc^), k — uj/c] nondipole effects in soft-x-ray photoemission are demonstrated 
via an experimental and theoretical study of angular distributions of neon valence photoelectrons 
in the 100-1200 eV photon-energy range. A newly derived theoretical expression for nondipolar 
angular distributions characterizes the second-order effects using four new parameters with primary 
contributions from pure-quadrupole and octupole-dipole interference terms. Independent-particle 
calculations of these parameters account for a significant portion of the existing discrepancy between 
experiment and theory for Ne 2p first-order nondipole parameters. 
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A mainstay of photoemission is the (electric-)dipole ap- 
proximation (DA), in which all higher-order multipoles 
are neglected Q. The range of validity of the DA has 
received renewed interest as recent experiments [^.^ un- 
covered breakdowns at progressively lower photon ener- 
gies. At high energies {hu! >5 keV), breakdown of the 
DA in photoionization is well-known, and a proper de- 
scription requires inclusion of many multipoles [Q. For 
soft-x-ray (Hlo <5 keV) photoionization, in contrast, first- 
order [0(fc)] corrections to the DA generally have been 
considered sufficient |5]. At these relatively low ener- 
gies, DA breakdown primarily leads to forward/backward 
asymmetries in photoelectron angular-distribution pat- 
terns. Especially striking have been observations of 
nondipole effects at energies below 1 keV a region 
in which the DA is usually considered valid. In the 
present work, an experimental and theoretical analysis 
of neon valence photoemission demonstrates a new, and 
unexpected, breakdown: significant second-order [0(fc^)] 
nondipole effects, primarily due to electric-octupole and 
pure-electric-quadrupole interactions, in low-energy pho- 
toemission. 

We begin with Cooper's 0{k) formula for the differen- 
tial photoionization cross section of a subshell (n, k) in 
a randomly oriented target using linearly polarized light 



+ {SuK + InK cos^ 6*) sin^ cos (/)|, (1) 

where cr„K is the photoionization cross section, /3„k de- 
scribes the angular distribution within the DA, and 
and 7„K are nondipole angular-distribution parameters 
characterizing the leading first-order corrections to the 
DA (mostly E2 — Ei terms). The angles 6 and </> are 
determined by the direction of the photoelectron relative 
to the photon-polarization e and photon-propagation k 
directions, respectively. The first two terms on the right 
of Eq. (P constitute the usual DA expression for the dif- 
ferential cross section, and the DA notion of a "magic 
angle" [9m = 54.7°, P2{cos9m) = 0] is preserved only in 
the 4) = 90° plane perpendicular to k. 

At this level of approximation, recent rare-gas ex- 
periments [^,^ observed significant modifications of 
photoelectron angular distributions from DA expec- 
tations, generally in good agreement with first-order 
independent-particle-approximation (IPA) calculations 
[^,0. The only exception is Ne 2p while measured 
values of 72s {52s is negligible when /32s — 2) agree fairly 
well with calculations, measured values of the combined 
parameter C2p (= 3J2p + 72p) are 30% larger than IPA 
predictions for energies near 1 keV. 

The same experiment also found /32p disagrees substan- 
tially with IPA calculations in this energy region, but is 
in close agreement with correlated calculations using the 



random-phase approximation (RPA) [^J^, thereby iden- 
tifying important electron-correlation effects well above 
the n — 2 thresholds This result led to speculation 
Q the discrepancy between measured and IPA-calculated 
C^2p values might also be due to interchannel-coupling ef- 
fects. However, subsequent first-order nondipole calcu- 
lations including electron correlation pT] disproved this 
notion; RPA values of Ne C,2p are in excellent agreement 
with the uncorrelated IPA results 

In this work, we explain much of this discrepancy be- 
tween theory and experiment for Ne Q2p- Beginning with 
theory, second-order [0(fc^)] corrections to the differen- 
tial cross section, which arise from interferences between 
El - E3, El - M2, E2 - E2, E2 - Ml, and Mi - Mi, and 
from retardation corrections to Ei — Ei amplitudes, are 
incorporated into Eq. (^i 
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r = ^{1 + {I3n. + A/3„«) P2(cos( 



+ (<5nK + 7iiK cos^ 0) sin cos 0+ r]„K P2{cos9) cos2(l) 

+ finK cos 20 -I- (1 + cos 2(/)) P4(cOS 6')| , (2) 

where the 0(A;^)-parameters A/3, 77, /i, and ^ are intro- 
duced |12|. Three of them satisfy the constraint rj + fj,+ 
^ = 0. Reference [|l2| contains complete formulae and a 
tabulation of first- and second-order parameters for all 
subshells of the rare gases helium to xenon. 

For experiment, we first present results for Ne 72s and 
C2p determined assuming only first-order corrections, em- 
bodied in i5 and 7, are needed to correctly interpret the 
data. The experiments were performed with an appara- 
tus designed to measure deviations from the DA ||l^ ; 72s 
and (^2p are determined using angle-resolved photoemis- 
sion intensities at 9m and both (p — 0° and </> — 90°. 
Figure 1 compiles old |3| and new values for C2p and 
72s (open squares) determined in this way. The solid 
curves represent 0{k) calculations which agree 

well with the 2s results, but disagree with the 2p results 
above 800 eV. 

To obtain 0{k^) predictions (dotted curves in Fig. 1), 
we first carried out numerical IPA studies of second-order 
corrections for neon. Wavefunctions for bound-state and 
continuum electrons were obtained from the radial Dirac 
equation in a modified Hartree potential. Values for all 
n = 2 angular-distribution parameters were calculated 
up to 2 kcV. Our results for 72s and (2p are in excel- 
lent agreement with previous nonrelativistic IPA calcu- 
lations Figure 2 shows values for A/3, 77, and ^ 

— —rj — /i) obtained from our second-order IPA calcu- 
lations. Primary contributions to these parameters come 
from El — E3 and E2 — E2 terms, with the octupole term 
contributing about 65% at 1 keV. A smaller contribution 
(w 10%) comes from the Ei — M2 term. 

For our measurement geometry jl^], Eq. (H) and our 
0{k^) calculations can be used to estimate the influences 
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of second-order effects on the analysis of experimental 
results for C,2p (= 3(52p + l2p)- Specifically, an effective 
value of C,, including these influences, can be defined as 



(C)cff 



(T(6'™,7r/2) 



7 + 3^ + V54(A^-7g/18) 
1-M 



(3) 



using the notation a{9,(f)) — ^{(),4>)- Using the results 
in Fig. 2, Coff for 2p and ^^s for 2s can be determined and 
compared to the results in Fig. f to identify influences of 
second-order nondipole effects. (We stress the measure- 
ments are fine, only the assumptions behind their analy- 
sis are suspect.) This exercise yields the dotted curves in 
Fig. 1, providing excellent agreement with 72s and clearly 
improved agreement with (^2p- The second-order correc- 
tions contained in ^cff account for much of the difference 
between the first-order theory and experiment for (^2p, 
demonstrating the first observation of 0{k^) effects in 
soft-x-ray photoemission. 

To confirm this unexpected finding, we made addi- 
tional measurements with our apparatus, which contains 
four electron analyzers in a chamber which can rotate 
about the photon beam. At a nominal angular position, 
two analyzers are at 9„i and — 0° in the plane perpen- 
dicular to the photon beam (0 = 90°), which we refer to 
as the dipole plane because first-order corrections vanish, 
while the other analyzers are positioned on the 35.3° cone 
in the forward direction with respect to the photon beam, 
with one of them at 9m and cj) = 0°. Photoemission in- 
tensities in the magic-angle analyzers are independent of 
P and can differ only because of nondipole effects. While 
the magic angle is no longer strictly valid when second- 
order effects are included, calculations show they can be 
unimportant in certain geometries (see below). 

New measurements were performed at ten rotational 
positions, yielding 20 angle-resolved intensities for Ne 2s 
and 2p photoemission at different angles 6 within the 
dipole plane, and 20 more at different angles 9 and 
around the nondipole cone. From the calculated results 
for A/32p in Fig. 2, direct second-order effects on f32p 
should be insignificant; A(32p ~ 0.005 near 1 keV, much 
smaller than our measurement uncertainties. Therefore, 
values of P2p determined from the dipole-plane spectra 
should agree well with DA calculations, i/ effects due to 
77, fj,, and ^ are negligible in the dipole plane. Using 
Fig. 2, the influence of these parameters on angle-resolved 
photoemission intensities can be predicted. In the dipole 
plane, we predict their effects will mostly cancel, and thus 
the excellent agreement []l0|| between experiment and the- 
ory for P2p is not surprising. 

In the nondipole cone, influences of the second-order 
parameters are superimposed on intensity variations due 
to the dipole /3 and the first-order 6 and 7 parameters. 



However, for both 2s and 2p, our calculations predict 
effects due to 77, fi, and ^ also mostly cancel around 
the nondipole cone. Furthermore, small residual effects 
around this cone are similar in sign and magnitude for 
2s and 2p, which is relevant because intensity ratios are 
the raw input for data analysis. Assuming no influence 
of second-order effects in the nondipole cone, we mod- 
eled the measured 2s/2p ratios around this cone using 
Eq. (|l|) to derive values for 72s and C2p- These results 
(solid circles in Fig. 1) agree extremely well with flrst- 
order calculations |]6|,[7 pT | , confirming our prediction of 
near cancellation of second-order effects in this geome- 
try. 

The above confirmation tests include two independent 
experimental methods to determine 72s and (2p- one re- 
lies on measurements at many angles in the nondipole 
cone, the other relies on comparison of magic-angle-only 
measurements in the dipole plane and the nondipole cone. 
For the former, second-order effects mostly vanish. For 
the latter, in contrast, the influences of 77, /i, and ^ on 
Ne 2p photoemission are expected to be opposite in sign 
for (f) = 0° and (j) = 90°, because of the cos(2(/)) terms in 
Eq. (^). Thus, second-order effects should be observable 
with the latter method, hence the differences in values 
for C2p derived using the two methods. 

As a demonstration of the influence of second-order 
nondipole effects on angle-resolved-photoemission inten- 
sities. Fig. 3 compares spectra taken with the two magic- 
angle analyzers. Figure 3a contains a neon photoemission 
spectrum taken at ?iijj=1200 eV, 9„i, and (f) = 90° in the 
dipole plane, where influences of S, and 7 vanish. In- 
cluded are fit curves showing modeled peak shapes and 
photoemission satellites to the left of the 2s peak. The 
overall fit (solid curve) matches the data very well, as 
indicated by the residual in Fig. 3b. 

This spectrum and fit are reproduced in Fig. 3c and 
compared to a nondipole-cone spectrum at 1200 eV, 9rm 
and = 0°. Intensity normalization between the spectra 
was achieved using 72s from Fig. 1, for which experiment 
and theory agree well. By inspection, the 2s/2p ratio is 
different in the two spectra. One possible explanation 
is a differential influence of first-order nondipole effects 
on the 2s and 2p intensities. As a quantitative test of 
this hypothesis, we derived the dotted region in Fig. 3c 
by multiplying the fit to the 2p peak in the dipole-plane 
spectrum by the expected differential effect on 2s and 2p 
peak intensities in the nondipole-cone spectrum deter- 
mined from IPA-/RP A- predicted values for 72s and C2p 
(see Eq. (|l|)). If first-order effects alone explain the ob- 
served variation of the 2s /2p ratio, then the dotted region 
should coincide with the 2p peak in the nondipole-cone 
spectrum. It does not, and thus the difference between 
the dotted region and the open-circle data (« 10%) is 
attributed to the influence of second-order effects. 

In conclusion, an experimental and theoretical study of 
valence photoemission from neon has demonstrated the 
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first observation of second-order (primarily Ei — and 
E2— E2) nondipole effects on photoelectron angular dis- 
tributions in the soft-x-ray region. A general expression 
for the differential photoionization cross section, includ- 
ing all contributions through second order, has been de- 
rived in a form convenient for comparison to experiment. 

This work was supported by NSF and DOE EPSCoR. 
AD and WRJ were supported in part by NSF grant PHY- 
99-70666. DWL acknowledges UNLV Sabbatical Leave 
support. The experiments were performed at the Ad- 
vanced Light Source at Lawrence Berkeley National Lab., 
supported by the DOE Materials Science Division, BES, 
OER under contract DE-AC03-76SF00098. 



* Present address: ITAMP, Harvard-Smithsonian Center 
for Astrophysics, Cambridge, MA 02138. 
[1] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of 
One- and Two-Electron Atoms ( Springer- Verlag, Berlin, 
1957). 

[2] B. Krassig, M. Jung, D. S. Gemmell, E. P. Kanter, T. Le- 

Brun, S. H. Southworth, and L. Young, Phys. Rev. Lett. 

75, 4736 (1995); M. Jung, B. Krassig, D. S. Gemmell, E. 

P. Kanter, T. LoBrun, S. H. Southworth, and L. Young, 

Phys. Rev. A 54, 2127 (1996). 
[3] O. Hemmers, P. Glans, D. L. Hansen, H. Wang, S. B. 

Whitfield, D. W. Lindle, R. Wehlitz, J. C. Levin, I. A. 

Sellin, R. C. C. Perera, E. W. B. Bias, H. S. Chakraborty, 

P. C. Deshmukh, and S. T. Manson, J. Phys. B 30, L727 

(1997). 

[4] H. K. Tseng, R. H. Pratt, S. Yu, and A. Ron, Phys. Rev. 

A 17, 1061 (1978). 
[5] D. W. Lindle and O. Hemmers, J. Electron Spectrosc. 

100, (1999) in press, and references therein. 
[6] J. W. Cooper, Phys. Rev. A 47, 1841 (1993); 42, 6942 

(1990). 

[7] A. Bechler and R. H. Pratt, Phys. Rev. A 42, 6400 

(1990); 39, 1774 (1989). 
[8] W. R. Johnson and C. D. Lin, Phys. Rev. A 20, 964 

(1979); W. R. Johnson and K. T. Cheng, ihid. 20, 978 

(1979). 

[9] M. Ya. Amusia and N. A. Cherepkov, Case Studies in 
Atomic Physics (North-Holland, Amsterdam, 1975), Vol. 
5, p. 155. 

[10] E. W. B. Bias, H. S. Chakraborty, P. C. Beshmukh, S. T. 
Manson, O. Hemmers, P. Glans, B. L. Hansen, H. Wang, 
S. B. Whitfield, B. W. Lindle, R. Wehlitz, J. C. Levin, I. 
A. Sellin, and R. C. C. Perera, Phys. Rev. Lett. 78, 4553 
(1997). 

[11] W. R. Johnson, A. Berevianko, K. T. Cheng, V. K. Bol- 
matov, and S. T. Manson, Phys. Rev. A 59, 3609 (1999). 

[12] A. Berevianko, W. R. Johnson, and K. T. Cheng, At. 
Bata Nucl. Bata Tables 74 (1999) in press. The parame- 
ters A, /tt, and v in this work have been replaced here by 
rj, fi, and ^ to avoid confusion with photon wavelength 



and frequency. 

[13] O. Hemmers, S. B. Whitfield, P. Glans, H. Wang, B. W. 
Lindle, R. Wehlitz, and I. A. Sellin, Rev. Sci. Instrum. 
69, 3809 (1998). 




500 



500 1000 
Photon Energy (eV) 
FIG. 1. Experimental and theoretical values of 72s and (^2p 
(3(52p + 72p) for neon. The effective quantities include sec- 
ond-order [0(A;^)] nondipole influences. See text for complete 
description. 
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FIG. 2. Nondipole parameters of O(fc^) for neon. 
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FIG. 3. Valence photoemission spectra of neon taken at 
1200 eV and Om = 54.7°. a) = 90° spectrum, including 
a fit. b) Residual of fit in (a), c) Same spectrum and fit as 
in (a) compared to a = 0° spectrum, d) Residual of fit to 
the = 0° spectrum (lower curve), and difference between 
the = 0° fit and the dotted region in (c) (upper curve). 
The hatched area is 2p photoemission intensity attributable 
to second-order corrections. See text for full explanation. 
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